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Recent measurements for F2(x,Q 2 ) have been analyzed in terms of the 'dynamical' and 'standard' parton 
model approach at NLO and NNLO of perturbative QCD. Having fixed the relevant NLO and NNLO parton 
distributions, the implications and predictions for the longitudinal structure function Fl(x,Q 2 ) are presented. It 
is shown that the previously noted extreme perturbative NNLO/NLO instability of Fl(x, Q 2 ) is an artifact of the 
commonly utilized 'standard' gluon distributions. In particular it is demonstrated that using the appropriate - 
dynamically generated - parton distributions at NLO and NNLO, Fl {x, Q 2 ) turns out to be perturbatively rather 
stable already for Q 2 > O (2 - 3 GeV 2 ). 
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1. Introduction 

The parton distributions of the nucleon are de- 
termined as a function of the Bjorken-a; variable 
at a specific low input scale Q = Qq mainly by ex- 
periment, only their evolution to any Q > Q be- 
ing predicted by QCD. In the standard framework 
(see, for example, |1I2| ) Qq is arbitrarily fixed at 
some value Qo > 1 GeV and the free parameters 
of the input distributions are varied iteratively 
without any constraint until the data and the pre- 
dictions yield a minimal \ 2 ■ I n this approach even 
negative gluon distributions in the small-a; region 
have been obtained [2], leading to negative cross 
sections like Fl(x,Q 2 ). Alternatively, within the 
dynamical parton model 3 4 5J, the distribution 
functions at Q > 1 GeV are QCD radiatively gen- 
erated from valence-like (positive) input distribu- 
tions at an optimally determined Qo = fJ. < 1 GeV 
(where 'valence-like' refers to a/ > for all in- 
put distributions xf(x,^ 2 ) ~ x a f (1 — x) b ' ). This 
more restrictive ansatz, as compared to the stan- 
dard approach, implies of course less uncertainties 
|4|5j concerning the behavior of the parton distri- 
butions in the small-a; region at Q > [i, which is 
entirely due to QCD dynamics at x < 10~ 2 . 

Following reference [Sj the perturbative stabil- 
ity of the structure function F L (x, Q 2 ) in the low 
Q 2 region, Q 2 < 5 GeV 2 , is studied within the 



framework of the dynamical parton model. For 
comparison the same analysis is repeated utilizing 
a set of 'standard' parton distributions previously 
determined in [7J. 

As pointed out in [2 8 9 10 the issue of the per- 
turbative stability of Fl(x, Q 2 ) in the very small- 
x region, x <, 10~ 3 , at the perturbatively relevant 
low values of Q 2 > (9(2 — 3 GeV 2 ), is important 
because it represents a sensitive test of the relia- 
bility of perturbative QCD. For the perturbative- 
order independent rather flat toy model parton 
distributions in [5j, assumed to be relevant at 
Q 2 ~ 2 GeV 2 , it was shown that next-to-next- 
to-leading order (NNLO) effects are quite dra- 
matic at j; < 10~ 3 (cf. Figure 4 of [8]). To 
some extent such an enhancement is related to 
the fact, as will be discussed in more detail in 
Section QJ that the third-order a 3 contributions 
to the longitudinal coefficient functions behave 
like xc L ~ —In a; at small x, as compared to 
the small and constant coefficient functions at 
LO and NLO, respectively. It was furthermore 
pointed out, however, that at higher values of Q 2 , 
say Q 2 ~ 30 GeV 2 , where the parton distribu- 
tions are expected to be steeper in the small-x 
region (cf. eq. (13) of 8J), the NNLO effects are 
reduced considerably. It is well known that dy- 
namically generated parton distributions [3] are 
quite steep in the very small-a: region already at 
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rather low Q 2 , and in fact steeper [1] than their 
common 'standard' non-dynamical counterparts. 
Within this latter standard approach, a full NLO 
(2-loop) and NNLO (3-loop) analysis moreover 
confirmed [219] the indications for a perturbative 
fixed-order instability observed in [8] in the low 
Q 2 region. Even additional resummations have 
been suggested [11) in order to remedy these in- 
stabilities, although such additional ad hoc small- 
x terms lack any quantitative theoretical basis 
within QCD. 

The plan of this contribution to the proceedings 
is as follows. Section 2 is devoted to a description 
of the formalism used for the NNLO and NLO 
analyses of deep inelastic scattering data within 
the framework of the dynamical parton model. In 
Section 3 the resulting dynamical parton distribu- 
tions of the nucleon and, in particular, their be- 
havior in the small-x region are discussed in com- 
parison with their 'standard' counterparts. The 
predictions for the longitudinal structure function 
are presented in Section [4] Summary and conclu- 
sions are given in Section 5. 

2. Theoretical framework 

The NNLO and NLO analyses presented here 
are performed in the common modified minimal 
subtraction (MS) factorization and renormaliza- 
tion scheme. Heavy quarks (c, b, t) are not con- 
sidered as massless partons within the nucleon, 
i.e. the number of active flavors appearing in the 
splitting functions and the corresponding Wilson 
coefficients is taken to be n/ = 3. This defines the 
so-called 'fixed flavor number scheme' (FFNS), 
which is fully predictive in the heavy quark sector: 
the heavy quark flavors are produced entirely per- 
turbatively from the initial light (u, d, s) quarks 
and gluons. It is nevertheless consistent and cor- 
rect to utilize the standard variable n t scheme for 
the function [T2] , 

In the MS factorization scheme the relevant 
structure function F2 as extracted from the DIS 
ep process can be, up to NNLO, written as 
[13|14|15] 

F 2 (x,Q 2 ) = 

F 2 + m (x, Q 2 ) + F 2iS (x, Q 2 ) + F%(x, Q 2 ,m 2 c ) (1) 



with the non-singlet contribution coming from 
the three active (light) flavors being given by 
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where <S> denotes the common convolution, 
a(Q 2 ) = a s (Q 2 )/A7r, C^ g (z) = 6(1 - z), C 2 NS 
is the common NLO coefficient function (see, for 
example, [16] ) and a convenient expression for the 

(2)4- 

relevant NNLO 2-loop Wilson coefficient C 2 ^ s 
can be found in [13]. The NNLO Revolution of 
the flavor non-singlet combinations = u + u — 
(d+d) = u v — d v and q£ = u+u+d+d— 2(s+s) = 
u v + d v + 4g — 4s, where q = u = d and s — s, is 
related to the 3-loop splitting function [17] P^j S , 
besides the usual LO (1-loop) and NLO (2-loop) 
ones, Pj^! and f^S > respectively [13118] . Notice 
that we do not consider sea breaking effects (u ^ 
d, s s) since the HERA data used, and thus 
our analysis, are not sensitive to such corrections. 
The flavor singlet contribution in (1) reads 



- F 2 , s (x,i 
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with£(a;,Q 2 ) = S g 



i(l+q) = u v +d v +4q+2s, 
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^2NS an< ^ ^ ne additional common NLO 
2g can be again 
Convenient expres- 



found in 



gluonic coefficient function C^ 1 ' 



for example. 



(2) (2) 

sions for the NNLO and C 2 ' g have been given 
in [14] and the relevant 3-loop splitting functions 
Pij \ required for the evolution of Y,(x,Q 2 ) and 
g(x,Q 2 ), have been derived in [19]. We have 
performed all Q 2 -evolutions in Mellin rt-momcnt 
space and used the QCD-PEGASUS program gU] 
for the NNLO evolutions, appropriately modified 
to account for the fixed n/ = 3 flavor number 
scheme with a running a s (Q 2 ). In NNLO the 
strong coupling evolves according to 

where /3 = 11-2//3, ft = 102-38//3 and ft = 
2857/2 - 5033//18 + 325/ 2 /54 and the running 
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a(Q 2 ) is appropriately matched at Q — m c and 
Q = Wfc. The values m c = 1.3 GeV and m& = 4.2 
GeV have been used, as implied by optimal fits 
[3] to recent deep inelastic c- and 6-production 
HERA data. 

The heavy flavor (dominantly charm) contribu- 
tion F% in (1) is taken as in [7|21j as given by the 
fixed-order NLO perturbation theory [22)23] . The 
small bottom contribution turns out to be neg- 
ligible. These contributions are gluon g(x,u. 2 F ) 
dominated where the factorization scale should 
preferably be chosen [24] to be /ip = 4m^ . It has 
been shown [4] that the resulting predictions are 
in perfect agreement with all available DIS data 
on heavy quark production and are furthermore 
perturbatively stable [24]. Even choosing a very 
large scale like \j? F = 4(Q 2 + 4m 2 ) leaves the NLO 
results essentially unchanged [25|26] . in particu- 
lar in the small-cc region. This stability renders 
attempts to resum supposedly 'large logarithms' 
(lnQ 2 /m|) in heavy quark production cross sec- 
tions superfluous. Since a NNLO calculation of 
heavy quark production is not yet available, we 
have again used the same NLO 0{a 2 s ) result. This 
is also common in the literature [27128129] and the 
error in the resulting parton distributions due to 
NNLO corrections to heavy quark production is 
expected [27] to be less than their experimental 
errors. 

3. Quantitative results 

For the present analysis the valence-like in- 
put distributions at Qo = fi < 1 GeV are 
parametrized according to [6|7|21] 

xq v (x, Ql) = N q7i x a - (1 - x) b "- (1 + Qj„ y/x 

+d Qv x + e qv x ' ) , 
xw(x, Qg) = N w x°" w {l — x) hl " (1 + c w \fx + d w x) , 

(5) 

for the valence q v = u Vl d v and sea w = q, g 
densities. Since the data sets utilized are insen- 
sitive to the specific choice of the strange quark 
distributions, we generate the strange densities 
entirely radiatively [3], starting from s{x,Q 2 ) ) = 
s(x,Qq) = 0, where Qo < 1 GeV. The normal- 
izations N Uv and Nd v are fixed by f Q u v dx = 2 



and J Q d v dx = 1, respectively, and N g is fixed 

via Jg x(E + g)dx = 1. The following data sets 
have been used: the small-x [3U] and large-x 
[3"T] HI Ff data; the fixed target BCDMS data 
[3J for Fl and F% using Q 2 > 20 GeV 2 and 
W 2 = Q 2 (± - 1) + m 2 > 10 GeV 2 cuts, and the 
proton and deuteron NMC data [33] for Q 2 > 4 
GeV 2 and W 2 > 10 GeV 2 . This amounts to a to- 
tal of 740 data points. The required overall nor- 
malization factors of the data turned out to be 
0.98 for HI and BCDMS, and 1.0 for NMC. We 
use here solely deep inelastic scattering data since 
we are mainly interested in the small-cc behavior 
of structure functions. The resulting parameters 
of the NLO and NNLO fits are summarized in 
Table [1] The corresponding dynamical gluon and 
sea distributions, evolved to some specific values 
of Q 2 > Q 2 , are very similar to the ones in [415] 
which were obtained from a global analysis in- 
cluding Tevatron Drell-Yan dimuon production 
and, at the NLO level, high-£r inclusive jet data 
as well. Furthermore, it should be mentioned that 
the NLO a s (M 2 ,) in Table Q] turns out be some- 
what smaller in fits based solely on deep inelas- 
tic structure function data [7|18|27|34|35] as com- 
pared to those which take into account additional 
hard scattering data [11214136] (for a recent sum- 
mary, see [37]). At NNLO the resulting a s (M§) is 
generally slightly smaller [37] (c.f. Table [T]) which 
is due to the fact that the higher the perturbative 
order the faster a s (Q 2 ) increases as Q 2 decreases. 

For comparison the results of a 'standard' fit to 
the same data [3] are presented in Tabled] where 
the gluon and sea input distributions in ([5]) do 
not vanish as x — > (a g .q <, 0) at Qq = 1.5 GeV 2 . 
Without loss of generality the strange sea at the 
input scale is taken to be s(x,Qq) = s(x, Qq) = 
0.5 q(x, Qq). The overall normalization factors 
of the data are 0.98 for BCDMS and 1.0 for 
HI and NMC, while slightly different values of 
the charm and bottom masses have been used, 
namely m c = 1.4 GeV and m& = 4.5 GeV. The 
standard NNLO and NLO fits are very similar 
to each other, corresponding to x 2 /dof = 0.989 
and 0.993, respectively, with the NNLO predic- 
tions for Fi falling slightly below the NLO ones at 
smaller values of Q 2 [7] . It should be emphasized 
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Table 1 

Parameter values of the dynamical NNLO and NLO QCD fits with the parameters of the input distribu- 
tions referring to |(SJ) at a common input scale Q 2 , = /i 2 = 0.5 GeV 2 optimal at both perturbative orders. 
Here y 2 wa s evaluated by adding in quadrature the statistical and systematic errors. 



NNLO NLO 





u v 


d v 
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Uy 


d v 


9 


9 


N 


0.6210 


0.1911 


0.4393 


20.281 


0.5312 


0.3055 


0.4810 


20.649 


a 


0.3326 


0.8678 


0.0741 


0.9737 


0.3161 


0.8688 


0.0506 


1.3942 


b 


2.7254 


4.7864 


12.624 


6.5186 


2.8205 


4.6906 


14.580 


11.884 


c 


-9.0590 


65.356 


2.2121 




-8.6815 


44.828 


-2.2622 


15.879 


d 


53.547 


1.6215 


7.7450 




54.994 


-5.3645 


21.650 




e 


-36.979 


-41.117 






-40.088 


-21.839 






X 2 /dof 






1.037 






1.073 










0.112 






0.113 







Table 2 

As TableQ]but for the standard NNLO and NLO QCD fits with the parameters of the input distributions 
referring to |5]) at a common input scale Qg = 1.5 GeV 2 . 



NNLO NLO 



U v 


d v 


9 


9 


u v 


d v 


q 


9 


N 0.2503 


3.6204 


0.1196 


2.1961 


0.4302 


0.3959 


0.0546 


2.3780 


a 0.2518 


0.9249 


-0.1490 


-0.0121 


0.2859 


0.5375 


0.2178 


-0.0121 


b 3.6287 


6.7111 


3.7281 


6.5144 


3.5503 


5.7967 


3.3107 


5.6392 


c 4.7636 


6.7231 


0.6210 


2.0917 


1.1120 


22.495 


5.3095 


0.8792 


d 24.180 


-24.238 


-1.1350 


-3.0894 


15.611 


-52.702 


•5.9049 


-1.7714 


e 9.0492 


30.106 






4.2409 


69.763 






X 2 /dof 

«a(M|) 




0.989 






0.993 








0.112 






0.114 







that the perturbatively stable QCD predictions, 
both in the dynamical and standard approaches, 
are in perfect agreement with all recent high- 
statistics measurements of the Q 2 -dependence of 
F2(x, Q 2 ) in the (very) small- x region. Therefore 
additional model assumptions concerning further 
resummations of subleading small- x logarithms 
(see, for example, [35]) are not required |17ll9j . 

The sea and gluon distributions resulting from 
both fits are shown in Figures Q] and [2] respec- 
tively. The dynamical NLO sea distribution has 
a rather similar small-a; dependence as the stan- 
dard one |4|7j : this is caused by the fact that 
the valence-like sea input in ([5|) vanishes very 
slowly as x — > (corresponding to a small value 
of a,q, dq ~ 0.05, according to Table [T]) and thus 
is similarly increasing with decreasing x down to 
x ~ 0.01 as the sea input obtained by a standard 



fit. On the other hand, the dynamically gener- 
ated NLO gluon is steeper as x — > than the 
gluon distributions obtained from the standard 
fits. Similar remarks hold when comparing dy- 
namical and standard distributions at NNLO. At 
NNLO the sea distribution xq is larger (steeper) 
than the NLO one, whereas the NNLO gluon dis- 
tribution xg is flatter as x decreases and, in gen- 
eral, falls below the NLO one in the small-x re- 
gion. It is evident from Figure that the NNLO 
gluon remains valencelike even at Q 2 — 2 — 4 
GeV 2 , i.e. decreases as x decreases; this is mainly 
caused by the dominant NNLO gluon-gluon split- 

(2) 

ting function Pg g which is negative and more 
singular as x — ► than the LO and NLO ones, 
P 9 (2) (x)^ilnl 0. 
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X 

Figure 1. The sea distributions xq{x,Q 2 ), where 
q = u = d, in the dynamical (dyn) and standard 
(std) parton model for two representative low val- 
ues of Q 2 . 

4. The longitudinal structure function 

F L (x,Q 2 ) 

In this section we turn to the perturbative pre- 
dictions for Fl(x,Q 2 ). Similarly to (Q]) one can 
write for Fl in the MS scheme, 

x~ l F L = C L , NS ®q NS 

+ ^ (C L , q ®qs + C L , g <8> g) + x- l F c L 

(6) 

where again ® in the nf = 3 light quark fla- 
vor sector denotes the convolution, q^s stands 
for the usual flavor non-singlet combination and 
qs = J2 q =u d sil + 9) ^ s the corresponding flavor- 
singlet quark distribution. We use the NLO ex- 
pression [22 23] for F£ also in NNLO due to our 
ignorance of the 0(al) NNLO heavy quark cor- 
rections. (Notice that F£ is a genuinely subdomi- 
nant NLO contribution to the total Fl, being less 
than 10% even at Q 2 = 2 GeV 2 and x = 1(T 5 




10" 5 10" 4 10" 3 10" 2 10" 1 10" 5 10" 4 10" 3 10" 2 10" 1 1 
X 



Figure 2. As in Figure [TJ but for the gluon dis- 
tribution xg(x,Q 2 ). 

and further decreases for increasing x. Further- 
more, the NNLO 3-loop corrections to F£ have 
been calculated recently [39] for Q 2 ^> m 2 , but 
this asymptotic result is neither applicable for our 
present investigation nor relevant for the major- 
ity of presently available data at lower values of 
Q 2 .) The perturbative expansion of the coeffi- 
cient functions can be written as 

C L 4a s ,x) = Y: (7) 

n=l ^ ' 

In LO, c« s = fx, c% s = 0, c« = 24x(l - x) 
and the singlet-quark coefficient function is de- 
composed into the non-singlet and a 'pure sin- 
glet' contribution, c^ q = c^ ns + c^ ps . Suf- 
ficiently accurate simplified expressions for the 
exact |40I4TI42| NLO and g3j NNLO coefficient 
functions Cj*\ and cf^\, respectively, have been 
given in [8]. It has been furthermore noted in 
[5] that especially for Cl g both the NLO and 
NNLO contributions are rather large over almost 
the entire x-range. Most striking, however, is the 
behavior of both C'l^ and Cl,q at very small val- 
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Figure 3. The individual light (u, d, s) quark con- 
tribution Fl to the total Fl in (J6j) in the dynam- 
ical (dyn) and standard (std) parton approach at 
NNLO and NLO for two representative low val- 
ues of Q 2 . The standard NLO results in the lower 
panel are similar for the CTEQ6 (anti)quark dis- 
tributions pp. Notice that, according to |6]), 



n 



n 



Ft- 



ues [8.44J of x: the vanishingly small LO parts 
(xCl\ ~ x 2 ) are negligible as compared to the 
(negative) constant NLO 2-loop terms, which in 
turn are completely overwhelmed by the positive 



NNLO 3-loop singular corrections xc 



O) 

L,i 



■In: 



This latter singular contribution might be indica- 
tive for the perturbative instability at NNLO [8], 
as discussed at the beginning, but it should be 
kept in mind that a small-a; information alone is 
insufficient for reliable estimates of the convolu- 
tions occurring in ([6|) when evaluating physical 
observables. 

We display the predictions for the convolutions 
of the individual light u, d, s quark {Fff) and gluon 



Figure 4. As in Figure [3] but for the gluonic con- 
tribution F[ to Fl in © with F[ = \x C L , g ® g- 



(F[) contributions in ([6]) in Figures [3] and [H re- 
spectively, at two characteristic low values of Q 2 . 
(Note that Fl + F[ = F L - F C L according to ©). 
Although the perturbative instability of the sub- 
dominant quark contribution in Figure [3] as ob- 
tained in a standard fit does not improve for the 
dynamical (sea) quark distributions, the instabil- 
ity disappears almost entirely for the dominant 
dynamical gluon contribution already at Q 2 ~ 2 
GeV 2 as shown in Figure 0] This implies that 
the dynamical predictions for the total Fl(x, Q 2 ) 
become perturbatively stable already at the rel- 
evant low values of Q 2 > 0(2 — 3 GeV 2 ) as ev- 
ident from Figure [5j in contrast to the standard 
results in Figure [5] In the latter case the stabil- 
ity has not been fully reached even at Q 2 = 5 



GeV where the NNLO result at x 



10" 



more than 20% larger than the NLO one. A simi- 
lar discrepancy prevails for the dynamical predic- 
tions in Figure [5] at Q 2 — 2 GeV 2 . This is, how- 
ever, not too surprising since Q 2 = 2 GeV 2 rep- 



Dynamical parton distribution functions 



7 



0.5 
0.4 
0.3 
0.2 
0.1 

0.4 
0.3 
0.2 
0.1 


0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 





: F L (x,Q 2 ) : 
Q 2 = 2 GeV 2 


, NNLO : 

: il >" nlo ; 

\ Q 2 = 3 GeV 2 : 




i i i i^'i 


\ Q 2 = 4 GeV 2 

; i i i i^'j 


\. Q 2 = 5 GeV 2 ; 

1 1 1 H^lMj 




\ Q 2 = 20 GeV 2 ■ 



10" 5 10" 4 10" 3 10" 2 10" 1 10" 5 10" 4 10" 3 10" 2 10" 1 1 
X 



0.5 

0.4 

0.3 

0.2 

0.1 



0.4 

0.3 

0.2 

0.1 


0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



: F L (x,Q 2 ) : 
\ Q 2 = 2 GeV 2 


, NNLO : 

: std nlo ; 

\ Q 2 = 3 GeV 2 ; 

1 1 1 "l^i 


; i i i i 


.\ </ = 4 GeV 2 


'"\ Q 2 = 5 GeV 2 ; 
1 1 1 "il^'i 


; i i i 


\ Q 2 = 20 GeV 2 ■ 



io -5 io" 4 io" 3 io" 2 io" 1 io" 5 io" 4 io" 3 io" 2 io" 1 i 



Figure 5. Dynamical parton model NNLO and 
NLO predictions for F L (x, Q 2 ) in @. 



Figure 6. As in Figure [5] but for the common 
standard parton distributions. 



resents somehow a borderline value for the lead- 
ing twist-2 contribution to become dominant at 
small x values. This is further corroborated by 
the observation that the dynamical NLO twist-2 
fit slightly undershoots the HERA data for F-i at 
Q 2 ~ 2 GeV 2 in the small-x region (cf. Figure 1 of 
[4]), which indicates that nonperturbative (higher 
twist) contributions to F2 become relevant for Q 2 
< 2 GeV 2 [31]. The NLO/NNLO instabilities im- 
plied by the standard fit results obtained in |2l9j 
at Q 2 < 5 GeV 2 are even more violent than the 
ones shown in Figure [6l This is mainly due to 
the negative longitudinal cross section (negative 
Fl(x,Q 2 )) encountered in [219] . The perturbative 
stability in any scenario becomes in general bet- 
ter the larger Q 2 , typically beyond 5 GeV 2 |8I2I9| . 
as shown in Figures [5] and [6] This is due to the 
fact that the Q 2 -evolutions eventually force any 
parton distribution to become sufficiently steep 
in x. 



For completeness we finally compare in Fig- 
ure [7] our dynamical (leading twist) NNLO and 
NLO predictions for Fl(x,Q 2 ) with a represen- 
tative selection of (partly preliminary) HERA- 
HI data [30I31I45I46] . Our results for F L , being 
gluon dominated in the small-x region, are in full 
agreement with present measurements which is 
in contrast to expectations |2|9j based on nega- 
tive parton distributions and structure functions 
at small values of x. To illustrate the manifest 
positive definitencss of our dynamically generated 
structure functions at Q 2 > /i 2 = 0.5 GeV 2 we 
show Fl{x, Q 2 ) in Figure [7] down to small values 
of Q 2 although leading twist-2 predictions need 
not necessarily be confronted with data below, 
say, 2 GeV 2 . As pointed out in [5], where also 
a study on the ±lcr uncertainty bands of Fl has 
been performed, future precision measurements 
of Fl could even distinguish between NLO results 
and NNLO effects in the very small-x region. 
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Figure 7. Our dynamical NNLO and NLO predic- 
tions for Fl at a fixed value of W = 276 GeV. The 
(partly preliminary) HI data [30 31 45 46] are at 
fixed W ~ 276 GeV. 

5. Summary and conclusions 

To summarize, recent deep inelastic data for 
the structure function , without the inclu- 
sion of any Fl data, have been analyzed in the 
dynamical and standard parton model approach 
at NLO and NNLO of perturbative QCD. In both 
approaches, perturbative QCD evolutions of par- 
ton distributions in the (very) small- a; region are 
fully compatible with all high-statistics measure- 
ments of the Q 2 -dependence of F2(x, Q 2 ) in that 
region. The results turned out to be perturba- 
tively stable, therefore additional model assump- 
tions concerning further resummations of sub- 
leading small- X logarithms are not required. 

Furthermore, the extracted parton distribu- 
tions have been used to predict Fl- It 
has been shown that the extreme perturbative 
NNLO/NLO instability of F L at low Q 2 , noted 
in [2l9ll0j . is an artifact of the commonly uti- 
lized 'standard' gluon distributions rather than 
an indication of a genuine problem of perturba- 
tive QCD. In fact it has been demonstrated that 
these extreme instabilities are reduced consider- 
ably already at Q 2 = 2 — 3 GeV 2 when utilizing 



the appropriate, dynamically generated, parton 
distributions at NLO and NNLO. It is interesting 
to notice, once again, the advantage of the dy- 
namical parton model approach to perturbative 
QCD. 
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